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[57] ABSTRACT 

A method of CW navigation is provided using four 
transmitters (TR1-TR4) operating at sufficiently low 
frequencies to assure essentially pure groundwave oper- 
ation. The transmitters are keyed to transmit constant 
bursts (I sec) in a time-multiplexed pattern with phase 
modulation of at least one for identification of the trans- 
mitters with the ability to identify the absolute phase of 
the modulated transmitter(s) and the ability to modulate 
low-rate data for transmission. The transmitters are 
optimally positioned to provide groundwave coverage 
over a service region of about 50x50 km for the fre- 
quencies selected in the range of 200-500 kHz, but their 
locations are not critical because of the beneficial effect 
of overdetermination of position of a receiver (R) made 
possible by the fourth transmitter. Four frequencies are 
used, at least two of which are selected to provide opti- 
mal lane resolution. All transmitters are synchronized 
to an average phase as received by a monitor receiver 
(M). 

8 Oaims, 5 Drawing Figures 
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LOW-FREQUENCY RADIO NAVIGATION 
SYSTEM 

ORIGIN OF INVENTION 5 

The invention described herein was made in the per- 
formance of work under the NASA contract and is 
subject to the provisions of Section 305 of the National 
Aeronautics and Sp^e Act of 1958, Public Law 85-568 
(72 Stat. 435; 42 use 2457). 

FIELD OF THE INVENTION 

This invention relates to radio navigation systems 
having broadcast communication capability, and more 
particularly to a CW radio navigation system using four ^ ^ 
displaced transmitters and one or more frequencies, 
with an accuracy of 5 meters over a service region of 
50X50 km. 

BACKGROUND OF THE INVENTION 20 

Low-frequency radio navigation systems have 
evolved from the early Loran (LOng RAnge Naviga- 
tion) system operating with pulsed frequencies between 
1800 and 2000 kHz and the Loran-C/D systems operat- 
ing at a frequency of 100 kHz. Each Loran station oper- 25 
ates in a pulsed sequential mode such that Only one out 
of three stations can be received at a time. During the 
transmission pulse period of the master station, both 
slave stations remain silent, after which the slave sta- 
tions will transmit pulsed energy in sequence. The time 30 
differences between pulses received from the slaves, as 
compared to the master, provide hyperbolic navigation 
in an acceptable manner for most purposes, particularly 
over a long range of about 1000 km. For navigation 
over longer ranges, the Omega system operating in the 35 
10.2 to 13.6 kHz band may be used, but accuracy is 
limited to 1 or 2 nautical miles at the present state of the 
art in predicting the phases of the Omega sky waves as 
functions of position and time. Omega transmits keyed 
(interrupted) CW signals from eight stations on several 40 
discrete frequencies utilizing multiplexing techniques. 

For navigation over very short ranges, such as to 100 
km from the transmitter, the pulsed format of Loran 
C/D is not necessary to discriminate the highly predict- 
able groundwave from the unpredictable skywave. 45 
Therefore, continuous wave (CVO signals, keyed in an 
Omega-like format, can be used. 

Instantaneous position and velocity may be deter- 
mined by phase comparison of CW transmissions from 
three or more synchronized transmitters. Two CW 50 
systems offering such capability are presently in opera- 
tion. One is the keyed CW Omega system referred to 
above, and the other is an unkeyed CW system known 
as Decca, which operates in the 70 and 130 kHz fre- 
quency band. A Decca receiver measures and integrates 55 
the relative phase difference between the signals from 
the transmitters which operate at different frequencies 
for identification. Tracor, Inc. of Austin, Tex,, 1ms stud- 
ied the use of a keyed CW system for river and harbor 
navigation, and Collins Telecommunications Systems 60 
Division of Rockwell International has amplified on 
that study. 

The present invention implements a keyed CW sys- 
tem to provide radio navigation for the Army’s pro- 
posed Mobile Automated Field Instrumentation System 65 
(MAFIS) using the highly predictable groundwave 
from a plurality of transmitters which transmits a single 
frequency or multiple frequencies in a time multiplexed 
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pattern. The choice of frequency is not critical, but 
frequencies in the range 200-500 kHz offer the widest 
groundwave coverage by minimizing skywave effects 
as a function of range. MAFIS is to be used for com- 
mand and control, communications, weapons-engage- 
ment, and navigation (porition and velocity) of up to 
2000 receivers, where each receiver inay be assigned to 
a soldier, tank or other mobile unit. This application 
requires 15-meter accuracy Over areas as large as 
50x50 km. Such a requirement cannot be met by any 
existing operational system. 

SUMMARY OF THE INVENTION 

In accordance with the present invention, four CW 
transmitters operating at a sufficiently low frequency to 
assure essentially pure groundwave operation are con- 
tinually keyed in a time-multiplexed pattern with bi- 
phase modulation (180® phase inversion), in such a pat- 
tern as to permit identification of the transmitters with 
the ability to identify the absolute (unmodulated) phase 
of their signals, and to suppress interference to the navi- 
gation from external sources. Four transmitters are used 
instead of the usual minimum of three to provide a more 
robust system in the sense of greater immunity from 
noise and geometrical dilution of precision. The time 
slot for each transmitter carrier burst is the same length. 
That simplifies phase measurements and phase tracking. 
Notwithstanding the phase modulation for transmitter 
identification, low-rate data can be transmitted in the 
form of short codes phase modulated on the carrier 
bursts synchronized with the bursts such that a receiver 
upon detection of data bits may infer the absolute signal 
phase so that positional tracking is unimpaired. 

The transmitters are optimally positioned to provide 
groundwave coverage over a service region of about 
50x50 km, so as to provide the maximum coverage 
with minimum skywave corruption of the groundwave 
for the frequency band selected. However, the locations 
of the transmitters are not critical because of the benefi- 
cial effect of overdetermination of position made possi- 
ble by the fourth transmitter. 

Although the four-transmitter system may operate 
with a time-multiplexed single frequency, it is prefera- 
ble to use four different frequencies time multiplexed 
from each transmitter for improving the overall signal- 
to-noise ratio of the system, and resolving ambiguity of 
position within the area of coverage. A phase control 
receiver monitors all four transmitters and continually 
corrects all to a common, equal value of phase, as re- 
ceived by the monitor receiver. The frequencies of the 
transmitters are so selected that they stand in ratios such 
as to provide optimal multiple-frequency lane resolu- 
tion. 

The novel features that are considered characteristic 
of this invention are set forth with particularity in the 
appended claims. The invention will best be understood 
from the following description when read in connection 
with the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a pictorial diagram of the transmitter net- 
work arrangement. 

FIG. 2 is a diagram of the time-division multiplex 
pattern of four frequencies transmitted by four transmit- 
ters in equal time slots, 

FIG. 3 is a diagram that illustrates lane ambiguity 
resolution using commensurate frequency signals. 
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FIG. 4 is a block diagram of an exemplary transmit- 
ter. 

FIG. 5 is a block diagram of an exemplary receiver. 

DESCRIPTION OF PREFERRED 
EMBODIMENTS 

Referring to FIG. 1 of the drawings, the preferred 
implementation of a radio navigation system utilizes 
four keyed CW transmitters TRl through TR4, each 
operating in the very low frequency band of approxi- 10 
mately 200 kHz to 500 kHz, for example near 300 kHz; 
a network phase-control monitoring receiver M with 
communications links to the transmitters to maintain 
phase synchronization among all transmitters; and any 
number of mobile navigation receivers R. In addition 15 
there may be a command control central station (not 
shown) and communications links to the monitor and 
transmitters for status monitoring, and for broadcast of 
low-rate data to the receivers via one or more of the 
transmitters. These various communications links (not 20 
shown) may be by cable or UHF radio. There may also 
be a radio communications link from each receiver to 
the command control central in order for each receiver 
to transmit its identification code and position data, 
either periodically or on demand from the command 25 
control central. 

As proposed for MAFIS, the navigation system is 
provided for monitoring the positions and velocities of 
the receivers. Therefore, the position and velocity data 
are not displayed at the receivers, but for other applica- 30 
tions, such as river and harbor navigation, such display 
could be provided on a chart overlay, or any other 
display, such as a printout of X and Y coordinates of the 
Universal Transverse Mercator (UTM) grid used by the 
Army and surveyors in general. 

The transmitters are optimally situated at locations 
approximately at the corners of a square having a 72 km 
diagonal for operation with essentially pure ground- 
wave coverage of a service region of 50 X 50 km. Actual 
transmitter locations may be based on minimizing per- 40 
turbing effects of vehicle traffic, terrain ruggedness, 
dielectric discontinuities such as deep fresh water, rock, 
etc., and reradiators, such as towers, power lines, and 
railroad tracks. While optimizing overall coverage, an 
attempt is made to obtain 15-meter accuracy in a prime- 45 
coverage subregion of the 50x 50 km service region. 
The choice of 300 kHz as the frequency, and 72 km as 
the diagonal of the service region is based upon maxi- 
mizing the extent of the service region for a specified 
amount of propagation error, principally the sky wave 50 
error. 

Propagation errors of the CW signals from the four 
transmitters pose the greatest difficulty in achieving this 
15-meter accuracy. The 285 kHz to 500 kHz range is 
selected to minimize sky wave contamination in the 55 
service region of 50x 50 km, and provide the greatest 
stability of propagation. The spacing of the transmitters 
is chosen on the basis of some analysis of the maximum 
range at which worst-case sky wave corruption of the 
groundwave would be negligible (40 dB below ground- 60 
wave). The maximum range obtainable under these 
considerations is approximately 72 km using frequencies 
in the 300 to 400 kHz subrange. 

The basic signal pattern is shown schematically in 
FIG. 2 for the time-multiplexed transmitters. The basic 65 
pattern has four time slots 1, 2, 3 and 4 and four frequen- 
cies fi, f 2 , fa and f 4 . In each time slot, a transmitter keys 
one of its frequencies on the air at constant amplitude, 
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and with constant phase. The transmitters transmit on 
these frequencies in cyclic time sequence- Transmitter 
TRl transmits on frequency fj in time slot 1, then on 
frequency f 2 in time slot 2, etc. Transmitters TR2, 3 and 
5 4 use the same pattern displaced by 1, 2 and 3 time slots, 
respectively. In each time slot, a transmitter keys one of 
its frequencies on the air at constant amplitude and with 
constant phase. The pattern repeats continuously. Note 
that every transmitter is transmitting a different fre- 
quency during each time slot so that all transmitters are 
on at all times for maximum power transmission into the 
service region for maximum total signal-to-noise ratio. 

Unlike the Omega system, all time slots are of equal 
length (0.25 sec), and so the frequency burst length is 
not to be used to identify the transmitters. Using equal 
length time slots greatly simplifies phase measurement 
and phase tracking. For transmitter identification, a 
different and unique technique is then provided in the 
following manner. One transmitter, such as transmitter 
TRl, is inverted in phase for one time slot every third 
basic pattern in the case of a single frequency embodi- 
ment in which a single frequency is time-multiplexed 
among all of the transmitters as shown for the first row 
of the pattern in FIG. 2. 

In the case of a multiple frequency system, more 
design freedom is available for structuring a transmitter 
identification system. For example, it would be possible 
to have all transmitters invert their phase in a predeter- 
mined time slot of every third basic pattern. Yet another 
convention might be to invert all phases of all frequen- 
cies during every third basic pattern. And to achieve a 
combined synchronization and interference suppression 
capability, the sequence of inversion over six patterns 

might be -h H — , where a minus designates phase 

35 inversion of all frequencies in a basic pattern. This bi- 
phase modulation sequence allows identification of the 
transmitters, identification of the absolute phase, and 
interference suppression. Shorter sequences of biphase 
modulation would serve to identify the transmitters, but 
that would not allow the receivers to infer the absolute 
phase (polarity) of the received signals. The result 
would be to effectively cut the tracking wavelength in 
half, thereby doubling the number of phase-modulus 
ambiguities. The interference suppression feature pro- 
vides the most powerful capability in the real world 
interference environment, and is therefore the preferred 
choice. 

A single-frequency embodiment using a frequency of 
about 300 kHz has a 1 km wavelength and 0.5 km lane 
width, which implies ambiguous position fixes spaced at 
approximately 0.5 km. Within the service region of 
50x50 km, the only available means for resolving these 
ambiguities are (1) a “consistency” method in which a 
search of the hyperbolic grid compares combinations of 
phases from four transmitters for geometrical consis- 
tency, and (2) using the leading edge slope of each 
narrowband carrier burst, as in Loran C/D. Lack of 
bandwidth on the transmitters may preclude the leading 
edge technique, and the “consistency” method may not 
be sufficiently reliable. 

Without some means for resolving lane ambiguity in 
a single frequency system, it would be necessary for the 
operator of a receiver to initialize the system to a 
known, or assumed, position within a half lane width. 
Ambiguity would then be resolved by selecting as the 
“correct” position the one that is closest to the assumed 
position out of all possible positions having the mea- 
sured phase angles of signals from all transmitters. In 
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the event there are two equally close possibilities, the 
ambiguity is resolved by some predetermined conven- 
tion, such as selecting the one to the north. Once that is 
done initially, the navigation problem is then just one of 
tracking. 5 

,In the four-frequency embodiment of the invention, 
however, lane ambiguity is resolved by the use of a 
commensurate relationship among the four frequencies 
such that 


mXi =iiX2=pX35i=^^4 


10 


20 


where fh, m p and q are iiiinimuih integers 6.^ » have no 
cbmmdn divisors) and Xi, X2, X3 and X4 are the Wave- 
lengths of the frequencies, prOVidM that the period of . . 
the commensuratene^ (mXi, nXi, pXa or qX4) is greater * 
than twice the Width of the service region^ 50 km in this 
example. However, for maximum robustness (iihmunity 
from noise and geometrical dilution of precison), the 
wavelengths (frequencies) should be so chosen that the 
commensurate “wavelength’^ of at least one pair of 
frequencies, e g. mXi with nXa, is the shortest one possi- 
ble to exceed twice the width of the service region. 

The measured phases of the separate commensurate 
frequency signals jointly contain information on the 
correct phase within the long commensurate “wave- 
length.” A mathematical likelihood function^ can be 
constructed and searched to find the correct phase. In a 
two-frequency embodiment, for example, ambiguous 
position in one of a large number of possible lanes can 
be resolved by effectively searching a predetermined 
table of phases of two commensurate frequencies fpr the 
closest match with a pair of phase measurements of the 
two commensurate frequency signals thmughout the 
entire commensurate “wavelength” of the two-fre- 
quency combination. The phase measurements of . the 
two measured commensurate frequency signals will 
have a closest match to the tabulated phases at the cor- 
rect lane. Lane resolution is thus achieved by using 
phase measurements from commensurate frequency 
signals and effectively looking up a table of predeter- 40 
mined phases for the argument of a pair which has the 
closest match to a measured pair. This tabular argun^ent 
provides lane resolution. How this is used for a position 
fix will be explained hereinafter with an example. If the 
two frequency signals selected are from the same trans- 45 
mitter, “radial” lane resolution is effected, and if from a 
pair of transmitters, “hyperbolic” lane resolution is 
effected. 

To construct a likelihood function for searching the 
table for the closest match, the differences between the 50 
measured phases and the tabular phases are computed. 
Then the sums of the squares of these differences are 
computed and a search is conducted for the minimum 
sum of the squares. In computing the differences be- 
tween the tabular phases and the measured phases, 55 
should any difference exceed 180*, its complement (360* 
minus the difference) is used in its stead. 

A two-frequency example will clarify how commen- 
surate frequencies aided by the likelihood function are 
used to resolve lane ambiguity. Assuming operation at 60 
near 300 kHz, frequency fi may be selected to be exactly 
300 kHz, and frequency fa may then be selected to be 
216/217 300 kHz. Both will then have lane widths of 
about one kilometer, i.e., the phase cycle for each re- 
peats itself about every kilometer so that phase measure- 65 
ments for determining position in the service area are 
ambiguous. This ambiguity can be resolved, however, 
by noting that the instantaneous phase of each relative 


to 'the other is unique over the cbmmehsuf ate wave- 
length consisting of 217 wavelengths at 300 kHz, i.e., 
throughout the entire span over which 217Xi coincides 
with 2I6X2, ' where Xi and X2 are the wavelengths at 
frequencies f] and f2. For every possible position in that 
span there Will be a different combination of phase mea- 
surements for the two commensurate frequencies. Con- 
sequently, having measured the phase of the signal at 

frequency fi and the phase <f>2 of the signal at frequency 
f2, there is only one lane for the signal of frequency or 
f2 Where the particular combination of <f»i and ^2 holds, 
as shown in an exaggerated form in FIG. 3 from Which 
it can be seen that for each point along the commensu- 
rate wavelength at 217Xi, there will be a unique pair of 
phase measurements of the two frequencies fi and f2. 
Predetermining and tabulating these phases With the 
accuracy of phase measurement that is possible in this 
navigation system (5 meters) will then permit resolution 
of lane ambiguity for the signal of frequency fi or fa. 
With this iane-resOlution capability built into the re- 
ceiver, no external input is required to obtain data from 
which a unique position fix in the 50X 50 km region can 
be obtained, and from which the phase tracking process 
can be initialized. 

To obtain a unique position fix, it may be noted that 
the lane resolution process determines the tabular argu- 
ment at which the tabulated phases most closely coin- 
cide with the measured phases. Thi^ tabular argument, 
so determined from the measured phases, is itself a 
“phase measurement” within the (long) commensurate 
“wavelength.” If, for example, “radial” lane resolution 
(described below) is done for two or more frequencies 
from each transmitter, the effect is mathematically in- 
distinguishable from having used a lower frequency for 
each transmitter. Thus, these “phase measurements” for 
three or more transmitters can be used to determine 
position. Because of the lower frequency mathemati- 
cally induced by lane resolution, the ambiguous position 
fixes (spaced approximately every half wave, for what- 
ever frequency is used) are more widely spaced. With a 
commensurate “wavelength” twice the span of the 
50 X 50 km region, unambiguous position fixing within 
this region is thereby obtained. 

The organization of the transmitters and receivers 
will now be described- FIG. 4 illustrates a functional 
block diagram of each transmitter which radiates the 
four keyed, time-multiplexed and phase-synchronized 
CW signals at commensurate frequencies fi through f4, 
A separate rubidium- or cesium-beam frequency stan- 
dard 10 is connected to each transmitter for open-loop 
synchronization of all transmitters. The monitor M 
(FIG. 1) continually corrects the phase of each trans- 
mitter to a common equal value of phase as received at 
the monitor M. It also synchronizes the keying of the 
transmitters in accordance with the pattern shown in 
FIG. 2. Use of a very stable frequency standard at each 
transmitter will minimize drift at the transmitters, and 
therefore minimize the rate at which the monitor M will 
need to transmit phase-control commands for maintain- 
ing phase synchronization among all transmitters. 

A binary counter and frequency synthesizer 12 is 
driven by the frequency standard 10 to generate a multi- 
ple n of the frequencies fi, f2, fa and f4. A suitable 
counter and frequency synthesizer is made commer- 
cially by Hewlett-Packard Co. (Model 6333). These 
frequencies are divided by n in binary counters 14 
through 17 to produce the desired frequencies fi 
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through f 4 . The multiple n is selected to be large, such 
as 2 ^^, in order that phase adjustments can be made in 
increments of 2^2 per frequency cycle by inhibiting or 
injecting counts in these counters synchronized with 
cycles of the frequency standard. These phase adjust- 5 
meats are made under command control from the moni- 
tor via a monitor receiver 18 and control unit 20. 

The monitor receiver 18 (coupled to the monitor M 
shown in FIG. 1 via cable or UHF radio) receives com- 
mands that are preceded by a transmitter identification 10 
code. Once the identification code is decoded, the moni- 
tor receiver 18 triggers the control unit 20 to advance or 
retard each of the frequency dividers 14 through 17 
separately a number of counts specified by the com- 
mand received from the monitor M. 

A local control loop slaves the antenna current to the 
desired reference phase in the frequency synthesizer 
using a transmitter current sensor 22 and control unit 
24, in a manner similar to that used in Loran C/D, 
which is to compare the antenna current phase at each 20 
of the frequencies to the phase of the frequency synthe- 
sizer. To facilitate implementation, a separate control 
unit may be used for each frequency, or a single, more 
complex unit may be time shared. In either case, closed 
loop control stabilizes the transmitter against local 25 
phase drift due to local conditions in the circuitry of the 
transmitter. 

The counter and frequency synthesizer 12 also drives 
a pattern format counter 26 which provides an output at 
the keying frequency to drive a time division multi- 30 
plexer 28, and a separate signal that is up (high) for 
every third period of the signal applied to the multi- 
plexer, for example, or for the -f H 1 — pattern 

referred to hereinbefore. The latter output is applied to 
a bank of biphase modulators 30, 31, 32 and 33 which 35 
selectively invert the phase of the signals at respective 
frequencies fi, f 2 , faand f 4 for identification of transmit- 
ter. This pattern format counter must also be synchro- 
nized with the pattern counter of all other transmitters. 
Once the transmitters are initially synchronized in this 40 
respect, the monitor M will maintain synchronization, 
just as it does phase synchronization of the frequencies 
transmitted. 

Also connected to at least one of the modulators, 
such as the modulator 33 for the frequency f 4 , is low- 45 
rate binary data to be biphase modulated on the carrier, 
although in practice, the low-rate data may be applied 
to all modulators simultaneously, as shown. In either 
case, the one transmitter biphase modulated for identifi- 
cation should not be used for this low-rate data commu- 50 
nication in order to permit receivers to continue track- 
ing the pattern and determine “sync” with the pattern as 
well as the absolute phase at all times. 

The biphase modulators are coupled to the multi- 
plexer 28 by waveshape filters 34, 35, 36 and 37 which 55 
smooth the digitally synthesized signals into sinusoidal 
signals at the fundamental frequencies fi, fz, f 3 and f 4 . 
This is in order to limit signal bandwidth. The output of 
the multiplexer drives a high-power, phase-stable ampli- 
fier 38 through an exciter 40, A multituned coupler 42 60 
couples the output of the amplifier to the antenna. The 
amplifier provides approximately 2 kW of power into a 
10% efficient vertically polarized antenna for radiation 
of about 200 watts. 

FIG, 5 is a functional block diagram of a four fre- 65 
quency navigation receiver of the open-loop type in 
which phase tracking and other navigational functions 
are accomplished by a programmed microprocessor. 
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Any signal from a low-frequency antenna 44 is picked 
off by an antenna coupler-preamplifier 46 (containing 
provision not shown for injection of a phase-shift cali- 
bration signal used to permit built in test as well as 
calibrate the phase shift as would be required for radial 
lane resolution). The preamplifier drives four parallel 
tuned receiver strips 48, 49, $0 and 51, each of which 
includes three sections A, B and C in cascade. 

The first section A is comprised of a filter-amplifier 
and limiter. The next section B includes a phase detec- 
tor and integrator. The last section C includes an ana- 
log-to-digital converter to provide digital phase mea- 
surements for processing. Phase detection and analog- 
to-digital conversion are synchronized by a timing gen- 
erator 52 having a temperature-compensated crystal 
oscillator 54 as its timing reference. The digital outputs 
are connected to a programmed microprocessor 56 via 
an input/output bus 58. Coupled to the bus are a display 
60 (if desired) and a coupler 62. The display is useful for 
calibration and diagnostics, but may also be used for 
display of position on a navigation overlay . The coupler 
is used to couple (to a cable or UHF radio) navigation 
data (coordinates of position X, Y; velocity X, Y; and 
other data, such as timing data) and to couple to a local 
printer, or the like, low-rate communication data re- 
ceived via the antenna 44 and decoded by the micro- 
processor 56. 

The coupler 62 is also used to couple in data, such as 
calibration data, e.g., propagation correction tables and 
coordinates of the transmitters and monitor at the time 
of setting up the navigation system. For that purpose, 
the coupler may include a keyboard, or other means to 
enter data. The microprocessor uses data from the strips 
48-51, and from the coupler 62, to generate position (X, 
Y) and velocity (X, Y) date according to any of the 
general and well known algorithms based on the 
method of least squares. 

In the following discussion for this algorithm, it will 
first be assumed that all of the transmitters are synchro- 
nized for equal phase at the transmitters for simplicity, 
following which the modifications to the algorithm for 
the method of synchronization of the transmitters at the 
monitor will be given. The following phase equations 
are shown in range units for simplicity: 

4>i=rXX y)+ 4 mofluloX (1) 

where r/ is the range as a function of X and Y coordi- 
nates from the ith transmitter plus some receiver offset 
time t measured when the receiver is referred to any of 
the transmitter locations, and where X is the wavelength 
at the frequency in question. A set of four equations 
given by equation (1) will be solved by a method of least 
squares. To implement the method of least squares, let a 
criterion, C, be defined by the following equation. 

4 (2) 

C= .2, n'Kha 10 + /]~4 >iF 
I 

The method of least squares seeks a value for variables 
X, Y and t such that the value of C is minimized. W is 
a weighting factor which, for the present invention is 
more nearly constant because the different keyed fre- 
quencies are for equal time slots. 

The method of least squares for the arrangement of 
synchronizing all transmitters to a common phase at the 
monitor, which is the preferred arrangement for the 
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present invention, will now be described. The measured 
phase <j>i at any constant, and at any location (X, Y), is 
related to the receiver’s position and instantaneous time 
offset (t/) with respect to the ith transmitter by the fol- 
lowing equation: 5 

4»/= mqduIoX (3) 

Note that this equation differs from equation (2) only in 
that time offset, t/, is different for the different transmit- 
ters because the monitor, located at known coordinates 
Xm, Ym so controls the transmitters that the phases of 
all transmitters are equal to a common, but unknown, 
phase value (with respect to the mobile receiver re- 
ferred to the monitor location) given by the following 
equation: 

Substituting t/from equation (4) into equation (3) yields 
the following equation: 

(5) 

Putting known quantities on the left (^/ is a measured 
quantity, and r/is the known range in wavelengths from 
the ith transmitter to the control monitor) equation (5) 
becomes: 

<f>i-hri(XM.rM)=ri(X,Y)+4* (6) 

Note that the right side of this equation differs from 
equation (1) only in that the unknown time offset t com- 
mon to all transmitters is replaced by an unknown phase 
common to all transmitters. The left side of this equa- 
tion may be viewed as a “modified measurement” of 35 
phase, 4>/. The measurement is modified only in that to 
each phase measured, the receiver must add the 
known quantity r, (XmiYm)- Once that modified mea- 
surement is entered, the data processing at the receiver 
for this method of synchronizing the transmitters to a 40 
common phase <j> at the monitor is the same as for the 
method of synchronizing all of the transmitters for 
equal phase at the transmitters. The least square solution 
minimizes 

4 (7) 

C ^ X Wi{[r,iX, y) + 4 .] - [4>| + riiXM Ym)]}^ 

/—I 

Note that the constant phase <f> in equation (7) to which 
all transmitters are synchronized at the monitor plays a 50 
role that is mathematically equivalent to the use of the 
time offset t in equations (1) and (2) for the first method 
of synchronizing the transmitters. The monitor could be 
placed at a precise point where the propagation time, t/, 
from each transmitter is the same, thereby simplifyingi#5 
the algorithm outlined in equations (3) to (7), but that 
would not be practical. Instead the monitor is optimally 
located to divide any error. 

Equations (1) to (7) effectively assume that the navi- 
gation receiver is stationary, i.e., the algebraic un- > 1 ^ 
knowns X, Y and t (or <j>) are constants to be obtained 
by solution of the equations. Should it be desired to 
accommodate a receiver moving at velocities X, Y or to 
accommodate a driftrate t (or ^) of the receiver’s time 
offset with respect to the transmitter system, the equa- 65 
tions (1) to (7) may readily be reformulated to include 
such velocities and driftrate. The method of least 
squares then produces algebraic solutions for X, Y, t, X, 
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Y, t in a manner similar to that used for stationary navi- 
gation. The generalization of equations (1) to (7) to 
accommodate a plurality of frequencies is equally 
straightforward. 

Once a receiver measures <J>/, the microprocessor 
effectively determines the space relationship of the re- 
ceiver to the transmitter and applies corrections using a 
three-dimensional position and temporal dependent 
propagation correction table, with interpolation be- 
tween tabulated values provided at approximately one 
kilometer intervals. Since the total phase delay is a 
function of the refractive index of the atmosphere, with 
a possible variation of as much as 20 to 30 meters of 
propagation over the total niaximum r^ge of ^m^ 

^ (70 km), it would be necessary to provide a plurality of 
tables for a plurality of refractive indices. These correc- 
tion tables are predetermined and stored in the receiv- 
ers. The command control central then periodically 
2 q issues an instruction to at least one of the transmitters to 
transmit to all receivers a data code indicating which 
table to use. 

The monitor M may be implemented with a receiver 
not unlike the field receivers. The receiver strips filter, 
25 amplify and detect the signal phases, and present to the 
microprocessor the phase measurements. As in the field 
receivers, the microprocessor identifies the transmitters 
from the biphase modulated pattern it receives. Having 
identified the transmitters, the microprocessor simply 
30 computes an average phase measurement, and then 
computes the corrections required to advance, or re- 
tard, the phase of the transmitters to that average phase. 
The corrections are then transmitted via cable or UHF 
radio to the respective transmitters. In the case of radio 
transmission, one command is transmitted for each 
transmitter on the same carrier with an operation code 
indicating phase correction, a transmitter identification 
code and the amount of correction for that transmitter. 
Each transmitter will then act on only the command 
having its identification code. The same technique 
would, of course, be used if all transmitters were con- 
nected to a common cable. 

What is claimed is: 

1. A method of low-frequency radio navigation com- 
prising the steps of 

transmitting from each of at least three spatially dis- 
placed transmitters at known locations a signal of 
at least one frequency selected from a range of 
about 200 kHz to about 500 kHz to assure operation 
with essentially pure ground wave, each transmit- 
ter being keyed on for an equal and constant burst 
length in a predetermined time-division multi- 
plexed pattern, and at least one transmitter being 
biphase modulated in a predetermined order of 
repetitions of said pattern to identify said one trans- 
mitter, and therefore all transmitters from the mul- 
tiplex pattern, 

phase synchronizing the transmitted signals from all 
transmitters to a common value of phase as re- 
ceived by a monitor at a known location notwith- 
standing the biphase modulation of one transmitter 
at predetermined intervals identification of all 
transmitters, thereby creating a region of zero 
error within the service area, 
receiving signals from said transmitters at an un- 
known position, and at said unknown position iden- 
tifying the absolute phase of said biphase modu- 
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lated signal in the multiplex pattern from the prede- 
termined order of the signals in the pattern, 
measuring the phase of each signal received at said 
unknown position from each transmitter, and 
from said phase measurements determining said un- 5 
known position by comparing measured phases 
with predetermined phs^es for all signals from all 
transmitters to determine the maximum likelihood 
of position. 

2. A method as defined in claim 1 wherein said moni- 10 
tor is positioned within an area bounded by said trans- 
mitters to divide any error due to variations in propaga- 
tion paths from the transmitters. 

3. A method as defined by claim 1 wherein the aver- 
age phase of signals received from all transmitters by 15 
said monitor is used as the common value of phase for 
phase synchronization of all transmitters. 

4. A method as defined in claim 1 wherein N transmit- 
ters are employed with as many frequencies as there are 
transmitters cyclically permuted by frequency and 20 
transmitter, where N is an integer greater than 2, 
whereby every transmitter is transmitting at all times on 

a different frequency, and all change frequency in uni- 
son in N time slots to produce a cyclic pattern of trans- 
mission from all transmitters. 25 

5. A method as defined in claim 1 wherein N transmit- 
ters are employed with up to N frequencies cyclically 
permuted by frequency and transmitter in a basic pat- 
tern and all transmitters are biphase modulated in a 
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predetermined order of repetitions of said patterns such 
that the number of inverted frequency patterns is equal 
to the number of uninverted frequency patterns for 
identification of transmitters and determination of abso- 
lute phase with interference suppression by allowing 
any interference to be integrated out of phase measure- 
ments. 

6. A method as defined by claim 1 wherein signal 
transmission from at least one transmitter is biphase 
modulated with low-rate data for communications of 
short messages and said transmitter biphase modulated 
for communication is other than said designated one 
modulated for transmitter identification. 

7. A method as defined by claim 1 wherein said signal 
transmission from said one transmitter biphase modu- 
lated for identification is modulated with low-rate data 
for communications of short messages during intervals 
between biphase modulation for transmitter identifica- 
tion. 

8. A method as defined by claim 1 wherein the signal 
transmission from at least one transmitter is biphase 
modulated with low-rate data for communication of 
short messages and the data code for communication is 
a synchronizable code with an equal number of biphase 
modulated periods representing encoded binary 1 and 0 
digits, thereby to allow interference from communica- 
tions to be integrated out of phase measurements for 
navigation. 

« « « * ♦ 
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